PTA/08-004 

A complete one-loop calculation of electroweak supersymmetric 
effects in ^-channel single top production at LHC 

M. Beccaria"'^, C. M. Carloni Calame'^, G. Macorini'^'^, 
E. Mirabella-^, F. Piccinini^, F. M. Renard'^ and C. Verzeernassi'^''^ 

° Dipartimento di Fisica, 
9P ' Universitd del Salento 

o ■ 



o 



(N 
O 
oo 
O 



X 



Via Arnesano, 1-73100 Lecce, Italy 



^ . ^ INFN, Sezione di Lecce, Italy 



INFN, Italy and School of Physics & Astronomy, 



University of Southampton, 

Oh| Southampton SOU IB J, England 

p ,■ , 

Q ■ Dipartimento di Fisica Teorica, 

^; 

Universitd di Trieste, 

cn 

' Strada Costiera I4, 

^ : 

0^ ■ Miramare (Trieste), Italy 



INFN, Sezione di Trieste 

f Max- Planck- Institut fur Physik 
Fdhringer Ring 6 D- 80805 Miinchen, Germany 
3 INFN, Sezione di Pavia, 

Via A. Bassi 6, 1-27100 Pavia, Italy 
^ Laboratoire de Physique Theorique et Astroparticules, 

UMR 5207 
Universitd Montpellier II, 
F-34095 Montpellier Cedex 5. 



1 



Abstract 

We have computed the complete one-loop electroweak effects in the MSSM for single top (and 
single antitop) production in the t-channel at hadron colliders, generalizing a previous analysis 
performed for the dominant dt final state and fully including QED effects. The results are quite 
similar for all processes. The overall Standard Model one-loop effect is small, of the few percent 
size. This is due to a compensation of weak and QED contributions that are of opposite sign. The 
genuine SUSY contribution is generally quite modest in the mSUGRA scenario. The experimental 
observables would therefore only practically depend, in this framework, on the CKM Wtb coupling. 

PACS numbers: 12.15.-y,12.15.Lk,13.75.Cs,14.80.Ly 
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I. INTRODUCTION 



The relevance of a precise measurement of single topproduction at hadron colliders was 
already stressed in several papers in the recent years [l]. A well known peculiarity of the 
process is actually the fact that it offers the unique possibility of a direct measurement of 
the Wtb coupling Vtb, thus allowing severe tests of the conventionally assumed properties of 
the CKM matrix in the Minimal Standard Model; for a very accurate review of the topics 
we defer to [2]. 

For the specific purpose of a "precise" determination of Vtb, two independent requests 
must be met. The first one is that of a correspondingly "precise" experimental measure- 
ment of the process. For the t-channel case on which we shall concentrate in this paper the 
CMS study [3] concludes that, with 10 fb^^ of integrated luminosity, one could be able to 
reduce the (mostly systematic) experimental uncertainty of the cross section below the ten 
percent level (worse uncertainties are expected for the two other processes, the s-channel 
and the associated Wt production, whose cross section is definitely smaller than that of the 
t-channel). 

The second request is that of a similarly accurate theoretical prediction of the observables 
of the process. In this respect, one must make the precise statement that, in order to cope 
with the goal of measuring Vtb at the few (five) percent level, a complete NLO calculation is 
requested. In the SM this has been done for the QCD component of the t-channel, resulting 
in a relatively small (few percent) effect jj]. The electroweak effects have been computed 
very recently at the complete one loop level within the MSSM for the dominant ub dt 
component of the process (to be defined in Section [IT]) 

One conclusion was that the genuine SUSY effect, for a set on mSUGRA benchmark points, 
was systematically modest, roughly at the one-two percent level. The SM contribution was 
computed in a preliminary way, i.e. only including the soft photon radiation to achieve 
cancellation of infrared singularities, ignoring the potentially relevant hard photon contribu- 
tion. Actually, the main purpose of 5| was that of investigating the possible existence and 
size of genuine supersymmetric effects, that would be essentially unaffected by the SM QED 
contribution. In this approximate approach, the one-loop SM contribution turned out to be 
sizable, of roughly ten percent on the total rate. This contribution should be considered as a 
perfectly known term, to be included in the theoretical expression of the rate and compared 
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with the corresponding experimental measurement to extract the precise value of the Wtb 
coupling Vtb- 

In fact, from the negative (for what concerns supersymmetric searches) result that the t- 
channel rate is not sensitive to genuine mSUGRA MSSM effects, adding the extra known 
feature that NLO QCD effects are small , of the order of five percent, and well under con- 
trol , one concludes that an extremely precise theoretical determination of the 
process would be possible, provided that a rigorous electroweak one-loop description were 
given. This requires two different steps: first, an extension of the calculation of sl to the 
seven remaining t-channel processes since the final dt state is only the numerically dominant 
one. Second, the additional calculation of the complete QED effect, including properly hard 
photon radiation. This is precisely the aim of this paper, whose goal will be that of offering 
a clean theoretical expression to be used for a significant measurement of Vtb- 
Technically speaking, this paper is organized as follows. In Section [TTl the definition of the 
eight considered processes will be given, and the main tasks that were fulfilled to perform a 
complete one loop electroweak calculation will be indicated. Since the problems to be solved 
were practically identical with those already met in p], our description will be whenever pos- 
sible quick and essential. In Section UTTt we shall expose the new calculation of the complete 
QED effects. In Section HVl we briefly discuss the one loop SUSY QCD effects that we have 
recomputed from scratch in the same scheme adopted for the electroweak corrections. In 
Section |V] we shall define the considered observables and show the results of our calculation, 
particular emphasis being given to the value of the total t-channel rate. A few conclusions 
will be drawn in the final Section I VII 



II. THE t-CHANNEL PROCESSES AT ONE ELECTROWEAK LOOP 

The complete description of the t-channel involves at partonic level four sub-processes 
for single t production: ub td, db tu, cb ts, sb —>■ tc, and the related four for the 
single i production. 

The starting point is the cross section for the ub — > td process with the complete set of 
one loop electroweak corrections (in the MSSM and SM). The 0{a) corrections to the 
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(unpolarized) differential cross section of this process read: 

spin 

where A^" is the tree level contribution to the amplitude of the partonic process uh dt 
while describes the EW one loop contribution to the same amplitude. The Mandelstam 
variables are defined as: 

t = {Pb-Ptf (2) 
u = {Pb-Pdf 

n 

The analytical expression for M.^ is available in literature (see for instance [5]). A^^ has 
been generated with the help of FeynArts |[9(], t^ algebraic reduction of the one loop 
integrals is performed with the help of FormCalc 



10| and the scalar one loop integrals are 



numerically evaluated using LoopTools [ll|]. We treat UV divergences using dimensional 
reduction while IR singularities are parametrized giving a small mass to the photon. 
The masses of the light quarks are used as regulators of the coUinear singularities and are 
set to zero elsewhere. 

UV divergences are cured renormalizing the parameters and the wavefunctions appearing 
in M.^ . In our case we have to renormalize the wavefunction of the external quarks, the 
mass of the W boson, the electroweak mixing angle and the electric charge. We use the on 



shell scheme decribed in Ref. 12|. This scheme uses the fine structure constant evaluated 
in the Thomson limit as input parameter. In order to avoid large logarithms arising from 
the runnning of a to the electroweak scale M^y, we slightly modify this scheme using as 
input parameter the Fermi constant Gp- We consistently change the definition of the 



renormalization constant of the fine structure constant following the guidelines of Ref. 13|. 



The unpolarized differential cross section for the process dh — ut can be obtained from that 
of the process uh —>■ td by crossing: 

d^^I^ut = ^ E 2^^{-^° *{s^u,u^ s)M\s ^u,u-.s)}. (3) 

spin 
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For the t production the cross sections can be calculated using the identities 



d<^I^tu = da2^ut (4) 

while the processes involving the second generation c, s and s, c quarks can be computed 
from the previous, simply replacing the masses of the external particles (and some masses 
in the loop corrections). 



III. QED EFFECTS 

In order to obtain physically meanigful observables one has to include the differential 
cross section for the process of t-channel single top production associated with the emission 
of a photon integrated over the whole photonic phase space. So we have to consider the 
partonic processes uh td'j, db tu'j, cb ts'j, sb — > tc'j, and the related four for the 
single t production. 

The unpolarized differential cross section of these processes has been obtained using two 
different procedures. In the first approach the amplitude has been generated and squared 
using FeynArts and FormCalc while in the second one the complete matrix element has 
been calculated with the help of FORM. The two methods are in mutual agreement. 
The phase space integration of the aforementioned differential cross section is singular in 
the region in which the photon is soft and in the region in which it is collinear to a massless 



quark. According to KLN theorem [ij] IR singularities and the collinear singularities related 
to the final state radiation cancel in sufficiently inclusive observables while the collinear 
singularities related to initial state radiation have to be absorbed in the redefinition of the 
Parton Distribution Functions (PDF). In order to regularize these divergences we use two 
different procedures: the Dipole Subtraction method and the Phase Space Slicing method. 
In the Subtraction approach one has to add and subtract to the squared amplitude an 
auxiliary function with the same asymptotic behaviour and such that it can be analytically 
integrated over the photon phase space. Among the different choices we use the function 



quoted in Ref. 



15|. In this Reference explicit expression for the subtraction 



for its analytical integration is obtained using the so called Dipole Formalism 



unction and 



3. 



The idea behind the Phase Space Slicing approach is to isolate the singular region of the 



6 



phase space introducing a cut off on the energy of the photon {AE = Ss\^/2) and on the 
angle between the photon and the massless quarks {6c)- In the regular region the phase 
space integration can be done numerically while in the singular region it can be performed 
analytically provided that the cutoffs are small enough. The form of the differential cross 
section in the singular region is universal and its explicit expression in the soft (collinear) 

n n 

region can be found in Ref. [1^ (|13|). 

As can be inferred from Fig. [2] the two methods are in good numerical agreement. 



IV. ONE LOOP SUSY QCD CORRECTIONS 



The one loop SUSY 



QCD corrections to t-channel single top production have been com- 



puted at LHC in Ref. [17|. We include these corrections re-computing them from scratch 
following the guidelines described in Sec. [Tll The only difference is that, following a stan- 
dard procedure in SUSY QCD, we treat UV divergences using dimensional regularization. 
Moreover in this case we have to renormalize only the wavefunctions of the squarks since 
the other renormalization constant do not have 0{as) corrections. These corrections are IR 
safe. 



V. OBSERVABLE QUANTITIES 

The differential hadronic cross section reads 



/I J r 

{da^i^ftis) + dal^.,,^{s) + rfaj- ,(.)) (5) 



(q. q')"^^0 

dL 



Where (g, g') = (u,d),{c,s),{d,u),{s,d). The differential luminosity has been defined ac- 



cording to Ref. [18| while tq = m'^JS and s = rS. c/aj (c/af^^°) are the 0{a^) (SUSY QCD) 
corrections to the differential cross section of the process X. 

As pointed out in Sec. Illlt initial state Collinear Singularities are not cancelled in the 
sum of virtual and real corrections and they are absorbed in the definition of the PDF. 
We use the MS factorization scheme at the scale = nit. Concerning the choice of the 
parton distributions set, we follow The calculation of the full 0{a) corrections to any 

hadronic observable must include QED effects in the DGLAP evolution equations. Such 



effects are taken into account in the MRST2004QED PDF [19[ which are NLO QCD. Since 
our computation is lea ding order QCD, we use the LO set CTEQ6L since the QED effects 



are known to be small 



20|. 



A. Numerical Results 

In this subsection we present our numerical results. We define Minv as the invariant mass 
of the (anti)top quark and of the light quark in the final state. Also, we denote by pt the 
transverse momentum of the (anti)top quark. We consider as physical observables the trans- 
verse momentum distribution the invariant mass distribution and two integrated 
observables derived from the previous: the integrated transverse momentum distribution, 
defined as the integral of from a minimum p™™ up to infinity: 

aipD = / -rrdp'^ (6) 

and the cumulative invariant mass distribution cr(M™^^), defined as the integral of the 
invariant mass distribution from the production threshold up to M™^^: 

r^^i-" da 

'J Tq inv 

For each observable we present the plots for the LO and NLO curves and for the percent- 
age effect of the NLO corrections to the observable [b = ^^^^^q^ x 100). Apart from the 
Standard Model results, we analyzed several supersymmetric mSUGRA benchmark points: 
in the figures we present the numerical results for the two representative ATLAS DC2 



mSUGRA benchmark points SUl and SU6 2l|]; in the other cases the results are similar. 



The SUl and SU6 input parameters generate a moderately light supersymmetric scenarios, 
where the masses of the supersymmetric particles are below the 1 TeV: the two physical 
spectra are quite similar, characterized by relatively heavy squark sector (562 — 631 GeV 
for the lightest squark in SUl and SU6 respectively), a light neutralino and a light chargino 
state. The typical masses for the supersymmetic Higgs are of order of 400 — 500 GeV . The 
main difference between the two points is the value of tan/5 input parameter: 10 for SUl 
and 50 for SU6. 
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1. SM Results 



In the Standard Model framework the behaviour of the four observables is shown in 
Figures El HI E] and [HI the left panels report the LO and NLO curves, and the right panels 
show the relative percentage effect of the one loop electroweak corrections: as one can see 
the curves for the NLO and LO are almost overlapping, and the global 0{a^) NLO effect is 
rather small. This is particularly evident considering the plot for the cumulative invariant 
mass distribution, where the percentage effect saturate to the ~ 1.5% for Mi„^ > 1500 GeV. 



2. MSSM Results 



The following figures show the analogous results for the transverse momentum distribu- 
tion the invariant mass distribution , the integrated pr ditribution (t(p™°) and the 
cumulative invariant mass distribution (j{M^^^) for the SUl (Fig. O [H [9l [TOl) and for the 
SU6 (Fig. [m [HI [131 [HD benchmark points. In the MSSM cases the NLO is defined as the 
sum of the electroweak part 0{a^) and of the SUSY QCD part; in each figure we show in 
the panel (a) the behaviour of the observable at LO and NLO; in (b) the global NLO effect 
(where "global" means 0{a^) plus SUSY QCD); in (c) and (d) the separate percentage 
contributions of the 0{a^) and SUSY QCD parts respectively. 

A first comment that can be drawn from the inspection of the plots is that the t-channel 
process is very weakly sensitive to the presence of the supersymmetric particle in the loops. 
The difference with the Standard Model case, below the percent level, is negligible for both 
the mSUGRA benchmark points, without appreciable differences between the two sets. 
The plots for percentage effect of the SUSY QCD corrections only (the (d) panel in Fig. 
[3 [HI [91 [lOjfor SUl and [IH [HI [El [H for SU6) show that the contribution of diagrams with 
virtual squark and gluinos is systematically small, below the 1% level, in agreement with 
Ref. 17|. The same conclusion holds for the pure electroweak supersymmetric contribution, 
being the 0{a^) effect due almost completely to the Standard Model part. 



3. PDF uncertainties 



An important source of theoretical uncertainty is given by the contribution of the para- 
metric errors associated with the parton densities. This could become particularly rele- 
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vant for single top channels, due to the presence of an initial state b quark, whose dis- 
tribution function is strictly related to the gluon distribution. We have studied the im- 
pact of such uncertainties on the transverse momentum distribution and on the in- 
variant mass distribution by using the PDF sets MRST2001E and CTEQ61E as 

. For each bin in the histograms the maximum and min- 
imum values are calculated starting from the central value according to the formula 



(j(central) ± 1/2^^^^20(15) (^(2^ ~ 1) ~ ^(2^))^, according to the prescription of Ref. |23| . 
The results are shown in Figure [151 The spread of the predictions obtained with the MRST 
set displays a relative deviation of about 2% or less, increasing on the large scale tails of 
the distributions (where the cross section is however very small), while the CTEQ set gives 
a larger uncertainty, of the order of 3-4%. This is due to different values of the tolerance 
parameter 2J], the latter being defined as the allowed maximum of the Ax^ variation w.r.t. 



the parameters of the best PDFs fit. Conservatively, we can associate to our predictions 
an uncertainty due to the present knowledge of parton densities of about 3%. It is also 
worth noting that the uncertainties obtained according to such a procedure are of purely 
experimental origin only {i.e. as due to the systematic and statistical errors of the data used 
in the global fit), leaving aside other sources of uncertainty of theoretical origin. 



VI. CONCLUSIONS 



We have computed in this paper the complete one-loop electroweak effect on various 
observables in the MSSM. The calculation has fully included QED effects. The overall 
result is that the one-loop effect is small, of a positive few percent in the total rate that 
we have considered as the first realistically measurable quantity. Technically speaking, this 
small number arises from a competition of negative weak contributions and positive QED 
terms. In the considered mSUGRA symmetry breaking scheme, the genuine SUSY effect in 
the considered benchmark points is systematically modest, at most of a one percent size. 
The values that we have obtained e.g. for the total rate should still be modified by the 
additional NLO QCD contribution. The latter is known, small and theoretically safe and 
could be easily added to our calculation. It will appear in a forthcoming paper that will 
provide an expression for the overall single top production at LHC, including the already 
existing calculations for the associated tW production (our paper and the QCD ones). We 
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have also given an estimate of the parametric errors associated with the present knowledge 
of the parton densities. The distributions studied in this paper are affected by a few percent 
PDF uncertainty. It is worth saying that this uncertainty is expected to be lowered once 
the LHC data will become available. In conclusion, a precise measurement of the t-channel 
rate appears as a perfect way to determine the value of the Vtb coupling, both within the 
SM and within the MSSM with mSUGRA symmetry breaking. We cannot exclude, though, 
that for different symmetry breaking mechanisms the genuine SUSY effect is more sizable. 
This question, that is beyond the purposes of this prehminary paper, remains open and, in 
our opinion, would deserve a special dedicated rigorous analysis. 
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FIG. 1: Born direct and crossed processes for single top production in the i-channel with first 
generation hght quark current. 
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FIG. 2: Lowest order partonic cross section for the process ub — > dt'j computed with the two 
different methods. The cuts for the Phase Space Shcing methos are 6s = Sc = 10~^. The quantity 
A is defined as A = a"'""'" 
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FIG. 3: Left Panel: We plot the LO (that is tree level) contribution and the NLO (that is tree level plus 
0{a^)) corrections to the transverse momentum distribution. 

Right Panel: We plot the percentage contribution of the 0{a^) corrections to the transverse momentum 
distribution; that is S = x 100. 

No cuts are imposed. Computation in the Standard Model framework 




FIG. 4: Left Panel: We plot the LO (that is tree level) contribution and the NLO (that is tree level plus 
0{a^)) corrections to the invariant mass distribution. 

Right Panel: We plot the percentage contribution of the 0{a^) corrections to the invariant mass distribution; 
that is (5 '^^NLo'^ X 100. 

No cuts are imposed. Computation in the Standard Model framework 



15 




FIG. 5: Left Panel: We plot the LO (that is tree level) contribution and the NLO that is tree level 
plus 0{a^) corrections to the integrated transverse momentum distribution a{p™™). We remind that this 
distribution is defined as the transverse momentum distribution integrated from a minimum p™'" up to 
infinity. 

Right Panel: We plot the percentage contribution of the 0{a^) correections to the integrated transverse 
momentum distribution; that is 5 = ^^nlo'^ ^ -'^OO- '^^ts are imposed. Computation in the Standard 
Model framework 
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FIG. 6: Left Panel: We plot the LO (that is tree level) contribution and the NLO (that is tree level plus 
0{a^)) corrections to the cumulative invariant mass distribution cr{Ml^^^). We remind that this distribution 
is defined as the invariant mass distribution integrated from threshold up to M™^^. 

Right Panel: We plot the percentage contribution of the 0{a^) corrections to the cumulative invariant mass 
distribution; that is (5 = ^^^lo'^' x 100- 

No cuts are imposed. Computation in the Standard Model framework 
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FIG. 7: (a) Wc plot the LO (that is tree level) contribution and the NLO (that is tree level plus 0{a^) 
plus SUSY QCD) corrections to the transverse momentum distribution. 

(b) We plot the percentage contribution of the 0{a^) plus SUSY QCD corrections to the transverse mo- 
mentum distribution; that is S — ^^nlq'^ ^ 100. 

(c) We plot the percentage contribution of the 0{a^) corrections to the transverse momentum distribution; 
that is (5 = X 100. 

(d) We plot the percentage contribution of the SUSY QCD corrections to the transverse momentum distri- 
bution; that is (5 = '^'^^J g*^-" X 100. 



No cuts are imposed. Computation in the SUl point 
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FIG. 8: (a) Wc plot the LO (that is tree level) contribution and the NLO (that is tree level plus C(ck^) 
plus SUSY QCD) corrections to the invariant mass distribution. 

(b) We plot the percentage contribution of the 0{a^) plus SUSY QCD corrections to the invariant mass 
distribution; that is S ^ ^^nlo'^ x 100. 

(c) We plot the percentage contribution of the 0{a^) corrections to the invariant mass distribution; that is 
<5=^xlOO. 

(d) We plot the percentage contribution of the SUSY QCD corrections to the invariant mass distribution; 
that is 5 = ^^^/^g^^ X 100. 



No cuts are imposed. Computation in the SUl point 



18 




FIG. 9: (a) We plot the LO (that is tree level) contribution and the NLO that is tree level plus 0{a^) 
plus SUSY QCD corrections to the integrated transverse momentum distribution cr(p™'"). (b) We plot the 
percentage contribution of the 0{a^) plus SUSY QCD corrections to the integrated transverse momentum 
distribution; that is S ^ ^^nlo" ^IQQ. 

(c) We plot the percentage contribution of the O(a^) corrections to the integrated transverse momentum 
distribution; that is (5 = x 100. 

(d) We plot the percentage contribution of the SUSY QCD corrections to the integrated transverse momen- 
tum distribution; that is 5 ^ ^"^nt n^^ ^ ^^O. 



No cuts are imposed. Computation in the SUl point 
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FIG. 10: (a) We plot the LO (that is tree level) contribution and the NLO (that is tree level plus 0{a^) 
plus SUSY QCD) corrections to the cumulative invariant mass distribution cr{Ml^^^). (b) We plot the 
percentage contribution of the 0{a'^) plus SUSY QCD corrections to the cumulative invariant mass distri- 
bution; that is <5 = '^^NLo'^ X 100. 

(c) We plot the percentage contribution of the 0{a^) corrections to the cumulative invariant mass distri- 
bution; that is S = X 100. 

(d) We plot the percentage contribution of the SUSY QCD corrections to the cumulative invariant mass 
distribution; that is 5 = ^'^'^^ x 100. 



No cuts are imposed. Computation in the SUl point 
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FIG. 11: (a) Wc plot the LO (that is tree level) contribution and the NLO (that is tree level plus 0{a^) 
plus SUSY QCD) corrections to the transverse momentum distribution. 

(b) We plot the percentage contribution of the 0{a^) plus SUSY QCD corrections to the transverse mo- 
mentum distribution; that is (5 ^^nlo'^ ^ 

(c) We plot the percentage contribution of the 0{a^) corrections to the transverse momentum distribution; 
that is (5 = X 100. 

(d) We plot the percentage contribution of the SUSY QCD corrections to the transverse momentum distri- 



bution; that is (5 = '^^^^g'^-" X 100. 



No cuts are imposed. Computation in the SU6 point 
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FIG. 12: (a) We plot the LO (that is tree level) contribution and the NLO (that is tree level plus C(ck^) 
plus SUSY QCD) corrections to the invariant mass distribution. 

(b) We plot the percentage contribution of the 0{a'^) plus SUSY QCD corrections to the invariant mass 
distribution; that is S ^ ^^nlo'^ x 100. 

(c) We plot the percentage contribution of the 0{a^) corrections to the invariant mass distribution; that is 
<5=^xlOO. 

(d) We plot the percentage contribution of the SUSY QCD corrections to the invariant mass distribution; 



that is 5 = ^^^/^g^^ X 100. 



No cuts are imposed. Computation in the SU6 point 
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FIG. 13: (a) Wc plot the LO (that is tree level) contribution and the NLO that is tree level plus 0{a^) 
plus SUSY QCD corrections to the integrated transverse momentum distribution (t(p™'"). (b) We plot the 
percentage contribution of the 0{a^) plus SUSY QCD corrections to the integrated transverse momentum 
distribution; that is S ^ ^^nlo'^ ^IQQ. 

(c) We plot the percentage contribution of the O(a^) corrections to the integrated transverse momentum 
distribution; that is (5 = x 100. 

(d) We plot the percentage contribution of the SUSY QCD corrections to the integrated transverse momen- 
tum distribution; that is 5 ^ ^"^nt n^^ ^ ^^O. 



No cuts are imposed. Computation in the SU6 point 
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FIG. 14: (a) We plot the LO (that is tree level) contribution and the NLO (that is tree level plus 0{a^) 
plus SUSY QCD) corrections to the cumulative invariant mass distribution cr{Ml^^^). (b) We plot the 
percentage contribution of the 0{a'^) plus SUSY QCD corrections to the cumulative invariant mass distri- 
bution; that is <5 = '^^NLo'^ X 100. 

(c) We plot the percentage contribution of the 0{a^) corrections to the cumulative invariant mass distri- 
bution; that is S = X 100. 

(d) We plot the percentage contribution of the SUSY QCD corrections to the cumulative invariant mass 
distribution; that is 5 = '^'^'^^g'^^ x 100. 



No cuts are imposed. Computation in the SU6 point 
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FIG. 15: Left Panel: PDF uncertainty (in percent) on the distribution (LO calculation). For each 
bin, the minimum and maximum deviations w.r.t the best fit PDF, as given by the MRST2001E set (solid 
lines) and by the CTEQ61E set (dashed lines), are shown. Right: the same as in the left panel, for j^f—- 
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